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The aim of the Pierre Auger Observatory is the investigation of the nature of cosmic ray 
particles at ultra-high energies. It can simultaneously observe the longitudinal air shower 
development in the atmosphere as well as particle densities on the ground. While there are no 
dedicated muon detectors, techniques have been developed to estimate the number of muons, 
A'^^, produced by air showers. Both, the longitudinal development, in particular the depth 
of the shower maximum, Xmax, and the muon content of air showers are highly sensitive 
to hadronic interactions at ultra-high energies. Currently, none of the available hadronic 
interaction models used for simulations of extensive air showers is able to consistently describe 
the observations of Xmax and made by the Pierre Auger Observatory. 



1 Introduction 

At the Pierre Auger Observatory^ extensive air showers induced by ultra- high energy cosmic ray 
particles with energies of 10^*^ — 10^'^ eV are investigated. The aim of the experiment is to solve 
the mysteries on the nature and the sources of these ultra-high energy cosmic ray particles. Since 
individual cosmic ray primaries and also ultra-high energy secondaries in the startup of the air 
shower cascade are interacting with our atmosphere at center-of-mass energies up to ~ 450 TeV, 
the Pierre Auger Observatory is also sensitive to the physics of hadronic interactions at center- 
of-mass energies far beyond the reach of the LHC. We demonstrate the current ability to test 
existing hadronic interaction models at these energies. 

The observatory consists of an array of 1600 water-Cherenkov particle detectors distributed 
on a surface of 3000 km^ combined with 24 fluorescence telescopes overlooking the atmosphere 
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Figure 1: Left panel: The total flux of cosmic ray particles'^'^ multiplied by the energy to the power of 2.5. 
Right panel: Measurements of (Xmax) compared to simulations^I. The uncertainty bands for proton and iron are 
indicating the spread caused by the hadronic interaction models QGSJetOI"^, QGSJetI]^, SlBYLL'^and Epo^. 



above the surface array in moonless nights. While the advantage of the surface array is a duty 
cycle of close to 100%, the fluorescence telescopes, which have a duty cycle of ~13%, provide 
an almost calorimetric measurement of the energy of the air showers as well as the position of 
the maximum energy deposition in the atmosphere, Xmax- 

The Pierre Auger Observatory combines the strengths of the two detection techniques, and 
can thus achieve an unpreceeded level of quality in data analysis. A good example is the 
measurement of the total flux, which is founded on the purely geometric acceptance of the 
surface array combined with the almost calorimetric energy reconstruction of the fluorescence 
telescopes^. This makes it the first cosmic ray experiment that is able to reconstruct the cosmic 
ray flux with an absolute minimum of systematic uncertainties induced by Monte Carlo. The 
resulting total cosmic ray flux is displayed in Figure [1] (left) , where it is shown together with the 
results of other experiments. The spectrum in the energy region of the Pierre Auger Observatory 
exhibits two important features. Around 10^®'^ eV there happens a hardening of the spectrum, 
and above ~ 10^^'^ eV a suppression is observed. The former is believed by most t heorists to 
be related to the transition from a galactic to an extragalactic origin of cosmic raysEIinilllll]. 
The latter is compatible with the GZK cutoff^^^IIll^ which is caused by interactions of cosmic 
ray particles with the extragalactic photon field, while it could also be related to the maximum 
acceleration energy of the cosmic ray source(s). 

On the other hand, for other analyses, as for example the reconstruction of the primary 
mass composition of cosmic rays, the prediction of air shower Monte Carlo simulations are a 
crucial ingredient The poorly constraint physics of hadronic interactions at the relevant 
ultra-high energies and phase space regions are currently preventing an unambiguous, or even 
self-consistent, analysis of the available data. At the same time it is demonstrated that air shower 
observables are very sensitive to hadronic interaction physics^^^. Thus, by using astrophysical 
constraints on the primary cosmic ray mass composition it is possible to study the features of 
hadronic interactions at ultra- high energies occuring in the startup of extensive air showers. The 
observables most sensitive to hadronic interaction physics are the depth of the shower maximum, 
^max; and the number of muons on the observation level, N^. 



2 The depth of the shower maximum 

The depth of the shower maximuni, -^max? is one of th.G most iiiiporta,nt observables of the 
telescope detectors. It is reconstructed by fitting a Gaisser-Hillas parameterization to the re- 
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Figure 2: Left panel: Measured muo nic signal, S^, in water-Cherenkov tanks obtained after subtraction of the 
electromagnetic signal component The lower markers are the predictions by QGSJetII for proton and the 
upper markers for iron induced air showers. The energy of the simulated showers was shifted up by 30 %. Right 
panel: Comparison of the muon reconstruction methods applied to Auger data. The muon number A^^ is defined 
as the relative muon number compared to the predictions of QGSJetII, and the energy scale ma y be varied 



within the systematic uncertainties, which are indicated as dashed vertical lines' 
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constructed energy deposit profiles-^. With the Pierre Auger Observatory a resolution better 
than 20g/cm^ for X^ax is achievedl^S]. in Fig. [T] (right) the measured mean Xmax versus energy 
is compared to predictions from air shower Monte Carlo simulations for different primary cos- 
mic ray particles and hadronic interaction models. The main features of the data of the Pierre 
Auger Observatory are a slight trend from heavier to lighter composition at low energies, while 
above lO^®'^ eV this trend is stopped, or maybe even reversed. The data can be described by all 
hadronic interaction models with a medium or mixed mass composition with A ^ 7. 



3 The number of muons at the observation level 



The Pierre Auger Observatory has no dedicated muon detectors. However, it is possible to 
indirectly infer the numb er of muons in the water-Cherenkov tanks of the surface array with 
several 

niethodstmHI]. 

The principle of one technique, which combines the data of the tele- 
scopes with that of the surface array, is the subtraction of the electromagnetic part of the air 
shower, as it is measured by the telescopes, from the signal in the surface array to yield the 
muon signal. The results of this method are shown in Fig. [2] (left). Even with the energy of 
the simulations scaled u p bv 30 %, which is just within the systematic uncertainties quoted by 
the Auger Collaboration^!^, the data are lying beyond the iron predictions. Similar results are 
obtained with all the other muon measurement methods on the data of the Pierre Auger Obser- 
vatory, which are summarized in Fig. [2] (right). The muon number larger than the predictions 
by iron primaries is by itself a strong indication of a deficient hadronic interaction modeling. 
Furthermore, comparing the results to the (ATmax) data exhibits a striking incompatibility of the 
interpretations with the existing interaction models. 

All the results discussed here are based on comparisons to simulations performed with 
QGSJetII. The choice of a different interaction model does not change the findings qualita- 
tively. Even the Epos model, which predicts the largest muon numbers'^, cannot fully account 
for the discrepancies. 



4 Summary 

It is shown that the two observables of the Pierre Auger Observatory that are most sensitive 
to hadronic interaction physics at ultra-high energies, Amax and A^, do not give a consistent 



picture when eompared to air shower Monte Carlo simulations. While the (Xmax) data is located 
well in between the predictions for proton and iron, the data is located outside of this interval 
and would favor a mass composition even heavier than iron. 

This result clearly points out the deficiency of existing hadronic interaction models to de- 
scribe cosmic ray data at ultra-high energies. The characteristics of hadronic interactions within 
the air shower cascade must be different from what is currently assumed in all the models. 

With improved constraints on the primary cosmic ray mass composition, observations of this 
kind can be used to study the relevant properties of hadronic interaction physics at center-of- 
mass energies up to ~ 450 TeV. 
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